Bose— Einstein condensate of kicked rotators with time-dependent interaction 
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A modification of the quantum kicked rotator is suggested with a time-dependent delta-kicked 
interaction parameter which can be realized by a pulsed turn-on of a Feshbach resonance. The mean 
kinetic energy increases exponentially with time in contrast to a merely diffusive or linear growth 
for the first few kicks for the quantum kicked rotator with a constant interaction parameter. A 
recursive relation is derived in a self-consistent random phase approximation which describes this 
superdiffusive growth of the kinetic energy and is compared with numerical simulations. Unlike in 
the case of the quantum rotator with constant interaction, a Lax pair is not found. In general the 
delta-kicked interaction is found to lead to strong chaotic behaviour. 
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The realization of Bose-Einstein condensation (BEC) 
of dilute gases has opened new opportunities for studying 
dynamical systems in the presence of many-body interac- 
tions. Previous investigations have addressed the effect 
of the nonlinearities due to the interactions on the dy- 
namical localization and on the integrability of the one- 
dimensional Gross-Pitaevskii equation (GP-equation) in 
the presence of a time-periodic delta-kicked external po- 
tential v(q,t) 00. As is well-known the classical coun- 
terpart of the quantum kicked rotator without interaction 
displays chaotic motion, leading to diffusive growth in the 
kinetic energy above a certain value of the kick-strength 
K of the external potential 0, . In the corresponding 
quantum kicked rotator after a few initial kicks, during 
which the mean kinetic energy increases roughly linearly, 
the motion becomes quasiperiodic (barring special cases 
where quantum resonances occur) with a dynamical lo- 
calization of the kinetic energy 0-0 • The infinite num- 
ber of conserved quantities, leading to the integrability 
of the quantum kicked rotator, are the probabilities with 
which the quasi-energy states are contained in a given 
initial state. The inclusion of the nonlinear GP-term to 
the Schrodinger equation in numerical simulations was 
found to give rise to chaotic behaviour and derealization 
of the kinetic energy 0,0- However, analytical consid- 
erations still show the existence of an infinite number of 
independent integrals of motion for the one-dimensional 
GP-equation with an arbitrary external potential • 
At the moment it seems difficult to reconcile this analyt- 
ical result with the results of the numerical simulations. 

In this paper we suggest and investigate a modification 
of the quantum kicked rotator by introducing a time- 
dependent interaction g(t) of delta-kicks. Furthermore, 
we conjecture that the numerically observed derealiza- 
tion of the quantum kicked rotator with constant inter- 
action parameter g could be a consequence of employ- 
ing discretizations of the nonlinear Schrodinger equation 
which do not preserve the property of integrability of the 



GP-equation with external potential in the continuum. 
The scaled GP-equation with external potential 
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V 2 q i/> + v(q,t) ^ + g H 2 V (1) 



is modified by introducing a time-dependent interaction 
parameter g(t) 



ihip = -^-V 2 .V ■ 
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v(q,t) $ + g(t) M 2 ^. 
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The time-dependent interaction can be achieved by a 
suitable time-dependent tuning to a Feshbach resonance 
[To|. By a short pulsed modulation of the magnetic 
field used for the tuning, it is possible to taylor a 
time-dependent coupling parameter of the form g(t) = 
?En=-oo 6(t — n) which we want to analyze here. Thus 
the BE-condensate, we consider in the following, receives 
additional delta kicks apart from that of the external po- 
tential already present for the quantum kicked rotator. 
In between two kicks, the BE-condensate, described by 
PJl follows free motion with the quantum kinetic energy, 
whereas in the original quantum kicked rotator with con- 
stant g (Eq. CQ), the motion is determined by the in- 
tegrable nonlinear Schrodinger equation (NLS) without 
an external potential. Therefore, numerical simulations 
become considerably more efficient and easier for Eq. (0) 
with delta-kicked interaction than for Eq. because 
the latter equation with permanent nonlinearity needs 
many integration steps between two kicks . 

The map for the wavefunction, stroboscopically taken 
after each kick t n , with delta interaction g(t) is given by 

</>(<?, = e-^ 9 ,W e ^/2-v^ (g)t + ) (3) 

V(q,t- +1 ) = K cos(q) + g \i>(q,t- +1 )\ 2 , (4) 

where the t+, t~ variables refer to times immediately 
after (+) and before (-) the occurrence of the n-th kick 
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Formally, Eq. J5J with a general interaction parameter 
g(t) can be transformed into Eq. by substitution of 

the wavefunction \/g{t)/ga ^>{q, t) = ip(q, t). In terms of 
the new wavefunction ip(q,t), the GP-equation acquires 
a spatially constant, but imaginary time-dependent po- 
tential part 



averaged kinetic energy (p 2 (t^ +1 )) can be related to the 
wavefunction of the previous kick 



-d q V{q,t n+1 ) ip(q,t+ +1 ) 



(16) 



th-iP = -—V 2 q ^ + v(q,t)ij + g \M 2 ^ (6) 
v(q,t) = v(q,t)+i^ (7) 

However, a Lax pair for the GP-equation with a complex 
potential v(q, t) and constant interaction go cannot be 
constructed as in @>[9] for a real external potential v(q, t): 
Using the following general ansatz with 2x2 matrices 
for the generator X(q, t) and T(q,t) of a Lax-pair 



X(q,t) = 
T(q,t) = 



-iki(q,t) y/goi>*(q,t) 
</g~o~<ip(q,t) ik 2 (q,t) 

( A(q,t) B(q,t) 
\C(q,t) D(q,t) 



(8) 
(9) 



with yet undetermined complex functions A, B, C, D 
and fci, fe, the compatibility of the pair of equations 



. dw dw 
n—— = X w n- 
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dq dt 
requires the following equations to be satisfied 
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C q + i(D - A) ip + 



thip = -i B q + i(D - A)iP* + 



'9o 
(fci 



(10) 



C (11) 
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/5o(C r-Bi;) (13) 
fg- Q (Cr -B$). (14) 

Comparing Eqs. (|11|) and l|12fl to (|6I7|I and its complex 
conjugate, the terms i(D— A)ip and i(D—A)ip* in I|11I12|) 
have to be identified with v(q,t) ip and v*(q,t) ip* and 
do not allow for an imaginary part of v(q, t) . Due to the 
missing of a Lax-pair, complete chaotic behaviour of the 
GP-equation with a general time-dependent interaction 
g(t) can be expected. This is confirmed by our numerical 
simulations. 

It is possible to derive an approximate recursion re- 
lation for the mean kinetic energy, just taken after kick 
time tZ 



(p 2 (4)) 



dq\d q ip(q,t+)f 



(15) 



and time (p 2 (tt+i))- The derivation is mainly 

based on the assumption that kick-to-kick correlations 
between ij){q, and ip(q,t+) can be neglected. Ap- 

plying the expression V(q,t~ +1 ) for the external poten- 
tial and interaction before the (n+l)-th kick (l.'il II) . the 
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d q V(q,t- +1 )r(q,4 +1 ) (17) 



x ( e -*/h-V( q ,t- +1 ) dq ^{q,t+) 



The last integral in 1|17[) contains only a single wavefunc- 
tion ip*(q, t^ + i) at time t n+ i so that, due to the assump- 
tion of independent phases between neighbouring kicks, 
this term approximately vanishes. In the second term of 
relation l|17|) . the spatial derivative d q ip(q,t+) is trans- 
formed by a unitary operator which cancels in the inte- 
gral so that the integrand becomes \d q tp(q, i+)| 2 , and in 



consequence the second term is (p 2 (t+)). Taking the spa- 
tial 
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tial derivative d q V(q, t n+1 ) in the first term of (p 2 (t^ +1 )} 



9 q V(q,t n+1 ) = -K sm(q) + g d q ip* (q,t n+1 ) ip(q,t n+1 ) 

+ jffct+il^W). ( 18 ) 

we only keep expressions with absolute values of wave- 
functions and their derivatives and obtain 



h 2 <^(C+i)> 
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+ 2g 2 J^dq\d q ^(q,t- +1 )\ 2 |<A(<Z, ^ +1 )| 2 |V(<?,i+ +1 )| 

1/(2tt) 1/(2tt) 

+ h 2 (p 2 (t+)), 

where the absolute values of ip(q,t~ +1 ) and Vfo^n+i) 
are replaced by the mean density l/(27r). Since the 
wavefunction ip(q,t~ +1 ) differs from ip(q,t+) only by the 
time development with exp{ih/2 ■ V 2 .}, the average of 
l 9 9^(?,^+i)| 2 can be replaced by \d q ip(q,t+)\ 2 . 

This results in a recursive relation between (p 2 (tn+i)) 
and (p 2 (t+)) 



(p 2 (4 + i)> = ^ + ( ' 



2tt 2 H 2 



(P 2 (4)). (20) 



In the case of vanishing interaction parameter g, (p 2 (t^)) 
increases linearly with K 2 /(2h 2 ), as expected for the 
initial classical diffusive regime of the original quantum 
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kicked rotator. However, the phenomenon of quantum lo- 
calization cannot be derived from the recursive relation 
which assumes no kick-to-kick correlation of the phases of 
the wavefunctions at time t n and t n +\ and neglects the 
correlations due to the presence of quasi-energy eigen- 
states. For finite interaction g quasi-energy states no 
longer exist and the assumption leading to (|2U[1 can be 
consistent. Indeed, due to Eq. I|2U[1 . the mean kinetic 
energy follows an exponential or superdiffusive growth, 
indicating strong chaotic behaviour, which makes the as- 
sumption of statistical independence of subsequent kicks 
self-consistent. The continuum limit of the recursion re- 
lation l|20(l also yields an exponential growth of (p 2 (t)} 



parameter g can be derived from a compatibility condi- 
tion with appropriate generators X and T for a Lax-pair 



{P 2 (t)} = — 5- 
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exp 
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where (p 2 (t — 0)) vanishes if the initial wavefunction 
if)(q,t = 0) = l/v2~7T is taken with zero momentum as 
in our simulations. Relation (|21|l for (p 2 (t)) interpolates 
between purely linear increase for small g and exponential 
growth for strong interactions. 

The comparison of the recursive relation for (p 2 (t^)) 
(I20|l and (p 2 (t)) <|21l) with numerical simulations yields 
a qualitative agreement on a logarithmic scale (see Figs. 
H to OJ. In Fig. we display on a logarithmic scale 
the mean kinetic energy for the delta-kick interaction 
g(t) = g X)n=^oo ^ — n ) w hh g = 5 and for the external 
potential v(q, t) = Kcos(q) X)"=-oo ~ n ) wu; h kick 
strength K = 1. Our estimate (J2UJ) and (|2TJ for (p 2 (t+)) 
gives a consistent value for the slope whereas the ab- 
solute values of {p 2 (t„)) differ by a constant factor from 
the numerical results. This deviation may result from the 
substitution of the average of the product of the densi- 
ties with the product of the average of the wavefunctions 
(cf. Eq. CHI)- The phenomenon of quantum localization 
as for vanishing interaction g = or small but constant 
g has not been obtained in our simulations with avail- 
able number of space points up to 2 16 on the periodic 
interval [0, 27r).The increase of the mean kinetic energy 
only ended at the maximum momentum squared which 
is limited by the spatial intervals. The complete missing 
of any sign of quantum localization and the superdiffu- 
sive increase of (p 2 (t)) are the main differences to the 
original quantum kicked rotator. In the case of constant 
interaction parameter g(t) — g (Eq. (Ql), a slow increase 
of the mean kinetic energy can be observed in numerical 
simulations 0,0- We conjecture that that this dereal- 
ization may be caused by a nonintegrable discretization. 
According to Ref. and our considerations, the obvious 
discrete version of the GP-equation Q is not integrable 
(ipn{t) = ip(q n ,t)) 



K 2 Tpn+i - 2V>» + ^ n _! 

' 2 (Ag)2 



+ V n (t) Ipn+g |-0"| 2 V'n 

and does not possess a Lax-pair, whereas the following 
discrete Schrodinger equation with constant interaction 



lhlp n 



h 2 ip n+ i - 2j> n + -0 ra _i 

■ 2 (Aq) 2 



(23) 



+ 2^ n (*) + Vn+l(t)) tpn + | \^n\ 2 (tpn+1 + V'n-l)- 

The mean kinetic energy derived from (|23|l and l|22|l can 
therefore be expected to behave quite differently, but lll'-'il) 
is, unfortunately, more difficult to simulate accurately 
than H22[). and has not been analyzed so far. 

Discrete versions of differential equations need not 
have in general the same integrability or chaoticity prop- 
erties as their continuous counterparts 9]. However, the 
modification of the quantum kicked rotator by a delta- 
kicked interaction, we have studied here, is free from such 
difficulties since it gives a unique map from the contin- 
uum equation J5J to the time development © of the 
wavefunction with only free motion between two kicks. 
The superdiffusive growth and strong derealization of 
the mean kinetic energy should be observable experimen- 
tally by a suitable tuning to a Fcshbach resonance. A 
realization of a ring structure for a BE-condensate, on 
which short potential kicks are applied, has been sug- 
gested |2j ■ Qualitatively, a stronger expansion rate of the 
condensate should be observed for the delta-kicked inter- 
action because of the exponential increase of the mean 
kinetic energy. 
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FIG. 1: Logarithm of the mean kinetic energy determined 
from the recursive relation 12U1 for vanishing initial mo- 
mentum of the wavefunction ip(q,to) = l/\/27r with a kick 
strength of the external potential of K = 1.0. \og{p 2 (t n )/2) 
is displayed for the first eight kicks t n versus the logarithm of 
the parameter g of the delta-kicked interaction, ranging from 
3 = 0.1 to g = 10.0. 
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FIG. 2: Results of the numerically obtained averaged kinetic 
energy on a logarithmic scale for the same parameters as in 
Fig- (Configurations of the axes are the same as in Fig. 
0) Qualitative agreement for the superdiffusive growth of 
{p 2 (t n )/2} is obtained with the exponential relation 1201 . 
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FIG. 3: Comparison of the exponential relation 112 H (dash- 
dotted upper line) for g = 5.0 of the delta-kicked interaction 
and the kick strength K = 1.0 with the corresponding nu- 
merical simulations (lower solid line). The slope of the ex- 
ponential increase of the kinetic energy, obtained by 1211 . is 
in good agreement with the computations whereas absolute 
values differ by a constant factor. 
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